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I based on ordered rooted trees, is studied. This Hopf algebra generalizes the 

Hopf algebraic structure of unordered rooted trees Ti-c, developed by Butcher 
in his study of Runge-Kutta methods and later rediscovered by Connes and 
Moscovici in the context of non-commutative geometry and by Kreimer where 
it is used to describe renormalization in quantum field theory. It is shown 
^ , that Ti.N is naturally obtained from a universal object in a category of non- 

' commutative derivations, and in particular, it forms a foundation for the 

study of numerical integrators based on non-commutative Lie group actions 
on a manifold. Recursive and non-recursive definitions of the coproduct and 
' the antipode are derived. It is also shown that the dual of Tipf is a Hopf al- 

\^ , gebra of Grossman and Larson. Ti. pf contains two well-known Hopf algebras 

' as special cases: The Hopf algebra He of Butcher-Connes-Kreimer is identi- 

I fied as a proper subalgebra of Hn using the image of a tree symmetrization 

" ■ operator. The Hopf algebra Tip oi the Free Associative Algebra is obtained 

from Ti.N by a quotient construction. 



O 

m 



Keywords Hopf algebra • ordered rooted trees • Lie group integrators 
Lie-Butcher series • Butcher group • Connes-Kreimer Hopf algebra • 
' Grossman-Larson Hopf algebra 



Mathematics Subject Classification (2000) 16W25 • 16W30 • 22E60 
37M99 • 65L05 • 65L06 • 81R60 • 81T15 



H. Z. Muntlie-Kaas 

Department of Mathematics, University of Bergen, N-5020 Norway 
E-mail: H.Munthe-Kaas@mi.uib.no 

W. M. Wright 

Department of Mathematics, La Trobe University, Victoria 3086, Australia 
E-mail: W.WrightQlatrobe. edu.au 



2 



1 Introduction 

Since Cayley [J] in 1857 and 100 years later Merson and followed shortly 
after by Butcher it has been known that rooted trees are extremely use- 
ful for structuring algebras of differential operators and elementary differen- 
tials. In 1972 Butcher ^ produced the far-reaching result that Runge-Kutta 
methods form a group. This was later named the Butcher group in the pa- 
per by Hairer and Wanner |17| . who made significant contributions to this 
theory. The Butcher group is defined on the dual of the tree space and it 
was pointed out by Diir that there exists a one-to-one correspondence 
between the Butcher group and the commutative graded Hopf algebra of 
unordered rooted trees. The Hopf algebra of unordered rooted trees has had 
far-reaching applications in various areas of mathematics and physics. In 1998 
a Hopf subalgebra was discovered by Connes and Moscovici fT^ during work 
on an index theorem in non-commutative geometry and by Kreimer 20.^ in 
the renormalization method of quantum field theory. Further collaborations 
between Connes and Kreimer [lOUll) have lead to other surprising results; 
notably a connection with the Riemann-Hilbert problem was established. 
Brouder [HE] realized that the mathematical structure of Connes-Kreimer 
was the same as that of Butcher. Grossman and Larson ^H] also developed a 
cocommutative graded Hopf algebra on a general class of rooted trees. It was 
shown by Foissy and Hoffman that the commutative Hopf algebra of 
Butcher and Connes-Kreimer was the dual of the cocommutative Hopf alge- 
bra of Grossman-Larson, which corrected the original result of Panaite |34| . 
Murua |3UII31j has developed series expansions of elementary differential op- 
erators and shown among other results that the logarithm of such a series is 
equivalent to the series expansions obtained from backward error analysis. 

Recently a great deal of interest has been focused on developing numerical 
methods which preserve geometric properties of the exact flow. In particular, 
Lie group integrators, which describe integrators that use Lie group actions 
on manifolds, were originally proposed by Crouch and Grossman |13| fol- 
lowed shortly after by Lewis and Simo |21II22| . Integrators of this type are 
now known as Lie group integrators; a survey of these methods is given in 
|19j . Series expansions for various classes of Lie group methods have been de- 
veloped; these expansions are generally used to analyze order. Munthe-Kaas 
|24II25 | constructed the order conditions for a special subclass of Lie group 
methods, where the computations are performed in a Lie algebra, which is 
a linear space. Later, [5^1 it was shown that the classical order conditions 
could be used along with a certain transformation. Owren and Marthinsen 
[55] developed the general order conditions for the Crouch-Grossman meth- 
ods with their analysis being based on ordered rooted trees. Recently, Owren 
|32j derived the order conditions for the commutator free Lie group methods 
[S], which were derived to overcome some of the problems associated with 
computing commutators. 

In this paper we aim to construct a commutative graded Hopf algebraic 
structure, which can be used to analyze the order of all Lie group meth- 
ods. The outline of this paper is as follows: In Section [21 we will introduce 
ordered trees and forests and describe some useful operations on them. We 
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will motivate the present Hopf algebra as a universal object in a category 
of non-commutative derivations, and also briefly discuss Lie-Butcher theory 
which will be treated in more detail in [55]. In SectionOwe develop the Hopf 
algebra of ordered trees, giving both recursive and non-recursive definitions 
of the coproduct and antipode, using certain cutting operations on ordered 
rooted trees. We show that the Hopf algebra described in this paper is the 
dual of the Grossman-Larson Hopf algebra, thus generalizing the result of 
Hoffman jl8| . Finally, in Section ^ we use a symmctcrization operator to 
provide an injective Hopf algebra homomorphism from the unordered trees 
into the ordered trees, establishing the former as a sub-algebra of the latter. 

2 Algebras of non-commutative derivations 

2.1 An algebra of trees 

In this section we will define an algebra N spanned by forests of ordered (and 
possibly colored) rooted trees. This algebra is a universal ('free') object in 
a general category of non-commutative derivation algebras, and plays a role 
in symbolic computing with Lie-Butcher series similar to the role of free Lie 
algebras 28,35 in symbolic computing with Lie algebras. 

Let OT denote the set of ordered colored rooted trees, and OF denote 
the (empty and nonempty) words over the alphabet OT, henceforth called 
the set of empty and non-empty forests. It should be noted that, unlike the 
classical Butcher theory, the ordering of the branches in the trees in OT is 
important, and likewise the ordering of the trees within the forest OF. 

The basic operations involved in building OT and OF are: 

— Create the empty forest I G OF. 

— Create a longer forest from shorter forests by concatenation, (a;i,Ci;2) >— > 



— Create a tree from a forest by adding a root node, : OF OT. In 
the instance where we wish to color the nodes using a set of colors X, we 
introduce an indexed family of root adding operations : OF — > OT 
for all i € T. The inverse operation whereby we create a forest from a tree 
by removing the root node is written B~ : OT OF. This operation 
extends to OF by B' {10002) = B' {loi)B- {L02) and ^"(1) = I. 

The total number of forests, with n nodes colored in i different ways, is 
defined by modifying the definition of the well known Catalan numbers 



See A000108 in "3^ for various combinatorial representations of the Catalan 
numbers. For a forest u) G OF we define the degree, #(w), as the number of 
trees in to as: 



L0lL02. 




#(1) = 0, 

#(S+(^)) = 1, 

#{lOiUJ2) =#{UJ1)+#{UJ2), 
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and the order, |w|, as the total number of nodes in all the trees of co as: 

|I|=0, 
\B+{uj)\ = l + U 
\coiu}2\ = |wi| + |a;2|- 

We let A'' = R(OF) denote the linear space of all finite M-linear combinations 
of elements in OF. This vector space is naturally equipped with an inner 
product such that all forests are orthogonal, 

(wi , a;2) = I J; = ' for all ui,U2 e OF. 

For a & N and lo G OF, we let a{ui) G M denote the coefficient of the forest 
CO, thus a can be written as a sum 

a = a{ui)w, 

where all but a finite number of terms are zero. The space of all infinite sums 
of this kind is denoted N* and is the dual space of N, that is 

N* = {a:N ^R: a linear}. 

We again let a{a) denote the value of a e N* on a G N. 

The operations , and concatenation extend to N by linearity and 
the distributive law of concatenation, that is 

B+(a;i + L02) = B+{loi) + B+ico^), 

U>{u)l + 0)2) = WWl + U)U>2- 

The vector space N with the concatenation product and the grading # forms 
a graded associative algebra N = ^j^^Nj, where Nj denotes the linear 
combination of forests with j trees. Alternatively, it is possible to grade this 
algebra using | • |; in this case, Nj denotes the linear combination of forests 
with the same number of nodes. 

Now we introduce a left grafting product which has the algebraic structure 
of a derivation. 

Definition 1 For a,a G N, define the left grafting a[a] € N hy the following 
recursion formulae, where t G OT and co,u! € OF: 

r[I] = 0, 

t[u)lo] = (t[w])w + c^{t[oj]), (la) 
r[B+H] = S+(rM) + i3+(ru;), 



l[a] = a, 

{Tu;)[a]=T[u;[a]]-{T[u;])[a]. ^ > 

The definition of left grafting is extended to the general case a[a] by bilin- 
earity. 
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It is useful to understand left grafting directly rather than via the recursive 
definition. From Hla|l we verify that if r S OT and u G OF then t[cij] is a sum 
of \uj\ words, each word obtained by attaching the root of t with an edge to 
the left side of a node of lu. 



• o o 

• ¥ i 



o 

i • o o 

i ¥ i 



o 

i up o 
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i 

loo 

• ¥ i 
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• i o 
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From ljlb|l we see that if ti,T2 G OT and u G OF then (Tir2)[a;] is obtained 
by first left-grafting T2 to all nodes of u; and then left-grafting ti to all the 
nodes of the resulting expression, except to the nodes coming originally from 

T2- 



O 
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Lemma 1 //ti, . . . , t/j G OT and lu G OF then (ti ■ ■ ■ Tk)[u!] is a sum of \uj\'' 
words obtained by, in the order j — k,k ~ 1, . . . ,1, attaching the root of the 
tree Tj with an edge to the left side of any node in lo. In particular we have 

oj[B+{l)] = B+{uj), for all uj G OF. (2) 

Equations (|la|l and ljlb|l imply that for any d & Ni and a,b £ N, we have 
the Leibniz rule and a composition rule of the form 

d[ab] = d[a]b + ad[b] (3) 

d[a[b]] = da[b] + d[a][b]. (4) 

Thus d acts as a first degree derivation on N. 

Definition 2 The Grossman-Larson (GL) product o : N(^$N — > is defined 

as 

B+ (wow) = [j[B+(w)], for aU G OF, 
and is extended to the general aod for a, a G A'^ by linearity. 

bmce LuoLu — B- (w[B+(5)]), the GL product can be understood by adding 
an invisible root to uj (turning it into a tree) , and left-grafting uj onto all nodes 
of B^{uj), including the invisible root. The root is then removed from each of 
the resulting trees, with the GL product resulting in a total of (|w| + !)#(") 
forests. Some examples of the grafting product from Definition ^ and the 
GL product from Definition [3 are given in Tabled In fa-ct, the theory of 
Grossman and Larson 16 is formulated on trees, not on forests of trees. To 
a forest lo in our terminology, they add a (proper) root to turn it into a tree. 
The definition of the GL product in (16) is modified accordingly. 

The GL product is an associative #-graded product, o : Nj®Nk Nj+k, 
satisfying for all a,b,c & N, aol = loa, and 

iaob)[c]^B-ia[B+ib)mc] 

^B-{B+{a[b]) + B+{ab))[c\ 
= {a[b] + {ab))[c] 
= a[b[c]]. 
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Table 1 The left grafting and Grossman-Larson products for all forests up to order 
four. The Grossman-Larson product is the dual of the coproduct in Tijv described 
in Section |3] 
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2.2 as a universal object 

Definition 3 Let D be an associative Z-graded algebra D = with 
associative product a,b i—>- ah, a unit I and grading ^{Dj) = j such that 
^{DjDk) — j + k. We call D a D-algebra if it is also equipped with a linear 
derivation (•)[■] : D®D — > D such that Q and Q hold for any d G _Di and 
any a,b G D. 

Define a D-algebra homomorphism as a linear degree preserving map J-^ 
between D-algebras satisfying for any a,b G D : 

T{ab) = T{a)T{h) (5) 

Ha[b]) = HaWm. (6) 

Proposition 1 Let N he the algebra of forests colored with a set X. For 
any D-algebra D and any map i <—>■ ft : T — > Di C D, there exists a unique 
homomorphism J- : N D such that J-{B^{T)) = fi. 

Proof From (01 and ® we find T{B+{uj)) = for any uj e OF. 

Together with (jsj and linearity, this shows that by recursion, we can extend 
J- to a uniquely defined homomorphism defined on all of N. □ 

This shows that iV is a universal object, free over the set 2, in the category 
of D-algebras. 

2.3 The algebra of 25 sections on a manifold 

As an example of a D-algebra, we consider an algebra related to the numerical 
Lie group integrators. Let g be a Lie algebra of vector fields on a manifold 
M and let exp : g T)iS{M), denote the flow operator. A basic assumption 
of numerical Lie group integrators [T??ll26ll33j is the existence of a g which is 
transitive (i.e. spans all tangent directions in any point on A^), and for which 
the exponential map can be computed efficiently and exactly. Transitivity 
implies that any vector field can be written in terms of a function f : M g. 
The goal of numerical Lie group integrators is to approximate the flow of a 
general differential equation 

- where f-.M^Q, (7) 

by composing exponentials of elements in q. The study of order conditions 
for Lie group integrators leads to a need for understanding the algebraic 
structure of non-commuting vector fields on A4, generated from /. 

Elements V G g are often called invariant or 'frozen' vector fields on 
Ai. These define first-degree invariant differential operators through the Lie 
derivative. Let V be any normed vector space and denote {Ai — > V) the set of 
all smooth functions from Ai to V, called the space of V-sections. For V £ g 
and ip £ {A4 V), the Lie derivative, ^^["0] G {A4 — + V), is defined as 



ip {exp{tV){p)) , for any point p E A4. 

t=o 
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For two elements V,W € q we iterate this definition and define the concate- 
nation VW as the second degree invariant differential operator = 
The linear space spanned by the 0-degree identity operator I[tp] = tp 
and all higher degree invariant derivations is called the universal enveloping 
algebra of g, denoted 25. This is a graded algebra with the concatenation 
product and degree #(I) = 0, #(0) = 1 and #{VW) = #{V) + #{W). 

Given a norm on the vector space 0, we consider the space of C5 sections"'^ 
{M ^<3). For two sections f,g e {M —>■ &) we define f[g] € {M ^ <d) 
pointwise from the Lie derivative as 

/[ff](p) = (/(p)[<7])(p), peM. 

Similarly, the concatenation on 65 is extended pointwise to a concatenation 

fg€{M^ 6) as 

{f9){p)=f{p)9{p), P&M. 
From these definitions we find: 

Lemma 2 Let f € {M ^ q) and g,h G {M ^ (&). Then: 

f[9h] = n9]h + gf[h], 
{fog)[h]=f[g[h]] = fg[h] + f[g][h]. 

Proof For p e M let V = fip) G Q. Then 



f[9h]{p) = 



d_ 
dt 
d_ 
dt 
_d 
di 



{gh){exp{tV){p)) 



t=o 



g{exp{tV){p))h{eMtV){p)) 



t=o 



g{eMtV){p))h(p) + g{p)h{eMtV){p)) 

t=o 

{f[9]h + gm){p), 



ifmmp) = 37 



{g{exp{tf){p))[h]{exp{tf){p)) 



t=o 



d 
dt 
d_ 
di 

if9)mp) + {f[9mh]{p). 



9{p)[h]{exp{tf){p)) + g{exp{tf){p))[h]{p) 



t=o 



□ 



Note the difference between fg and fog. In the concatenation the value 
of g is frozen to g{p) before the differentiation with / is done, whereas in the 
latter case the spatial variation of g is seen by the differentiation using /. 

^ Thus 25 is a trivial vector bundle over A4, g a trivial sub-bundle and the tangent 
bundle TM is a non-trivial sub-bundle of fl. 
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LemmaElshows that {Ai ©) is a I?-algebra. Thus if we, for every i G X, 
pick a vector field fi G {A4 g) then there exists a unique homomorphism 
: N ^ {M^(8) such that T{B+{I)) = The images of the trees J^{t), 
for T G OT, are cahed the elementary differentials in Butcher's theory (see 
[nj) and the images of the forests ^(w), for cu G OF, are cahed elementary 
differential operators in Merson's theory (see 12iij ). 

2.4 Elements of Lie-Butcher theory 

To motivate the algebraic structures of the next section, we briefly introduce 
some elements of Lie-Butcher theory. This theory is the non-commutative 
generalization of the classical Butcher theory and is the general foundation 
behind the construction of order conditions for Lie group integrators. Various 
aspects of this theory have been developed in [^24II25II27II33| . A comprehen- 
sive treatment is given in |29| . 

With the vector space N being defined as the finite R-linear combinations 
of OF, we now let N* denote the space of infinite R-linear combinations (se- 
quences), or the algebraic dual space of N. All the operations of the previous 
paragraphs extend to N* by local finiteness, see [HHl. For example, the inner- 
product in 12.1|l is extended to a dual pairing (•, •) : N* xN -^R, where the 
computation (a, b) = J2uj£OF <^{^)b{^) always finite, since b is required to 
be finite. 

Consider the homomorphism J-' introduced in Section 12.31 extended to 
a homomorphism of infinite series J-: N* {A4 ©), where {A4 25) 
should now be understood as a space of formal series. The series might not 
converge, but all definitions make sense termwise, and any finite truncation 
yields a proper ©-section. In classical (commutative) Butcher theory the 
image of a G N* is called an S-series, see Murua ISHj- Similarly, we define an 
LS-series as an infinite formal series in {Ai — » 0), given by ^ 

LS(a) = J2 ^'"''a('^)->^(w). (8) 

Classical Lie series on manifolds is a generalization of Taylor series, where the 
fundamental result is the following 'pull-back formula': Let / G {A4 ^ g), be 
a vector field and exp(/) : Ai A4,he the t = 1-fiow. For any g G {A4 g) 
we have (see PP) that 

oo 

g(exp(/)(p)) = J2 -]J'i9]{p) ^ Exp(/)[g](p), 

where = I and f^g] - /[• ■ • ' ' ■ ] = (/« ' • ' °/)[.9]- Note that if 
/ = J-(^ ), then f^ = J-{^ o • ■ ■ o» ), thus the operator exponential Exp(/) = 
Ejlo /Vj! is a LS-series. 

Two special cases of LS series are of particular importance: A LS-series 
LS(a) is called logarithmic or algebra-like if LS(a) G (A^ — > g) represents 



^ See comments at the end of Section |1] on the chosen normalization. 
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a vector field, and it is is called exponential or group-like if LS(a) is the 
(formal) operator exponential of a logarithmic series. A logarithmic LS-series 
is the non-commutative generalization of a B-series, named the Lie-Butcher 
series ^ . 

Note that if ti,T2 £ OT, then a = T1T2 — T2T1 is a logarithmic series since 
it represents the commutator of two vector fields. More generally, a series a S 
N* is logarithmic, if and only if, all its finite components belong to the free Lie 
algebra generated by OT. A Hall basis for this space is characterized in |27| . 
Reutenauer presents several alternative characterizations of logarithmic 
and exponential series. 

We find the characterization in terms of shujfle products particularly use- 
ful. The shuffle product U : N^N — > iV is defined for two forests as the 
summation over all permutations of the trees in the forests while preserving 
the ordering of the trees in each of the initial forests, and is extended to 
N(SiN by linearity. It can also be recursively defined in the asymmetric way 
lU Lo = Lu Ul = Lu for any forest uj G OF, and if uji = tiVi and 0^2 — T2V2 for 
Ti,r2 G OT and Vi,V2 G OF, then 

(tiVi) U {T2V2) = Ti{vi U (T2t;2)) + T2((tiUi) U V2). 

The shuffle product is associative and commutative, for all LUi,uj2,io^ G OF 
we have 

(wi U UJ2) Li LU3 = LUi Li {l02 Li [J3), 

UJl Li LU2 = LiJ2 Li UJl. 

Table El gives some simple, but nontrivial, examples of the shuffle product. 
Lemma 3 l'-{5f A series a G N* is logarithmic if and only if 
a{l) = 

a{uJiULU2)^0 for all UJ I, UJ2 € OF 
A series a G N* is exponential if and only if 
a{l) = 1 

a{uji Ua;2) = a(tJi)a(a;2) for allLUi,LU2 G OF. 

The LS-series of an exponential series a G N* represents pull-backs, 
or finite motions on A4. They form a group under composition with the 
Grossman-Larson product, which is the generalization of the Butcher group 
to the case of non-commutative group actions. Consider the GL product as a 
linear operator o : N*®N* — » N* defined by o[u!i®u!2) = u!iou!2- To compute 
the composition ao(3 for general a,(3 G N*, it is useful to introduce the dual 
of o, the linear map An : N N^N to be defined in Section mi Using 
Corollary n we find 



(ao/3)(tj) = (o(a®/3),w) = (a®/3,Z\jv(t^)) = ^ a(wi)/3(w2)- 
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Table 2 All nontrivial examples of the shuffle product for all forests up to and 
including order four. The shuffle product is the dual of the coproduct in the 
Grossman-Larson Hopf algebra. 

As an illustrative example, we read from Tabled that 

(ao/3)(. I ) = a(. I )m + M. . )+«(• ) + "(. )/?(. . ) + I )■ 

The inverse in the group is found from the antipode, a linear map Sn '■ 
N N defined in Section O It can be shown that the dual of the antipode 
5^ : N* N* also defines the inverse in the group: 

a~^{uj) = a(5Ar(w)) = (a, S'Ar(w)} = {S'^{a),u), 
for all exponential a e N* . 

3 Hopf algebras 

In this section we will study a commutative graded Hopf algebra TCn oi 
ordered trees. The coproduct Z\jv in the algebra is defined by recursion for- 
mulae, and at a later stage we will show that this An is the dual of the GL 
product thereby establishing the connection between this Hopf algebra and 
the algebra of the Butcher group. 

We begin by briefly reviewing the definition of a Hopf algebra, see |27] for 
details. A real associative algebra ^ is a real vector space with an associative 
product /z : A(s^A — > A and a unit u : R — > ^ such that /i(a®u(l)) = 
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^{u{l)(S)a) = a for all a G A. The dual of an algebra is called a coalgebra, C, 
which is a vector space equipped with a coassociativc coproduct A : C ^ C®C 
and counit e : C — > M. A bialgebra S is a linear space which is both an algebra 
and also a coalgebra structure such that the coproduct and the counit are 
compatible with the product, in the sense that 

e(/u(a;i(8)W2)) = ;u(e(wi)(8)e(a;2)), (9) 
Z\(M(wi(8)a;2)) = (/U(g)M)(/(g)T(g)/)(zi(a;i)(8>Zi(a;2)), (10) 

where T{coi(E)U!2) = W2<8iwi is the twist map. Let End(B) denote all linear 
maps from B to itself. We define the convolution ★ : End(B) — > End(S) as 

iA-kB){a) = tx{{A(SB)A{a)), for A, B e End(B) and a G S. (11) 

Let / e End(S) denote the identity matrix. An antipode is a linear map 
S e End(B), which is the two-sided inverse of the identity matrix under 
convolution, with the antipode satisfying 

{I*S){a) = {S*I){a)=u{e{a)), for all a e B. (12) 



Definition 4 A Hopf algebra W is a bialgebra equipped with an antipode. 



3.1 The Hopf algebra of ordered trees 

We will study a particular Hopf algebra based on the vector space of ordered 
forests A'' = IR(OF), where the coproduct is defined by the following recursion. 

Definition 5 Let An : N — > N(SiN be defined by linearity and the recursion 

An{I)= lOl, 

^ _ (13) 

+ Zijv(w) U •(/®B+)Z\Ar(S-(T)), 

where r = Bf{u)) S OT and a;,a; G OF. The linear operation 

U- : N®N®N®N — »■ N®N is a shuffie on the left and concatenation on the 

right, satisfying 

(wl^Ti) U •(W2®T2) = (Wi U Wi)(8)(tiT2). 

Note that letting w = 1 yields the special recursion formula for a tree r: 
An{t) = T®1+ {I®B+)An{B- [t)). 
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Theorem 1 Let Ti.N be the vector space N = M(OF) with the operations 

product: /iAr(a(g)6) = a U 6, (Shuffle product) 

coproduct: An, (Definition\^ 
unit: MAr(l) — I, 

coumt: e^H-|o^ ^i^^ 

Then Ti^ is a Hopf algebra with an antipode Sn given by the recursion 
§Jv(II) =1, 

(14) 

where r — [uj) G OT and lo,uj ^ OF. In particular 

%n{t) = -m ((§jv«'/)(/«>S+)Z\jv(5~(t))) . 

Proof The coassociativity of the coalgebra wih be estabUshed once we have 
shown that the algebra and coalgebra are compatible. From the fact that 
eAr(ti;) = for all ui G OF\{I} and that the shuffle product of two scalars is 
just standard multiplication, we immediately get l^. To show (|1U|I we find 
it convenient to introduce the linear operation LU : N®N(^N(^N N(^N 
with the shuffle product both on the left and the right, satisfying 

The compatibility condition (|10|) is now equivalent to 

An{i^iUuj2) ^ An {uJi)UJ An {i02)- (15) 
To simplify the notation we use the fact that 

An{cJt) — LJT®1 + Z\Ar(wT). 

Let uji — TiVi and lu2 — T2V2 for ti, T2 S OT and vi,V2 G OF, now using the 
recursive definition of the shuffle product and substituting the expression for 
the coproduct of an ordered forest yields 

An[viTi U V2T2) = An{{vi U T2V2)ti) + A{{tiVi U U2)T2) 
= {AN{v^)VnAN(v2T2)) U ■ZAr(ri) + (zAw (uin) L±J Z\Ar (uz)) U ^n{t2) 
+ {viTi U V2)t2®1 + [vi U V2T2)ti®1 

= {An{vi) LijZAr(u2r2)) U ^n{ti) + {An{viTi) UJ'An{v2)) U ^n{t2) 
+ {AN{vi)inV2T2®l) U •ZAr(Tl) + {An {V2) VIA VlTl®l) U ^n{t2) 
+ {viTl U U2T2)(8lI 

= ViTi U V2T2®1+^n{viTx) LU {V2T2®1) + ~An{v2T2) LU {viTi®!) 
+ An{viTi) LU An{v2T2) 
= An{viTi) LU A{V2T2)- 
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We now have the tools needed to show that the coalgebra is coassociative, 
which foUows from the fact that the coproduct A]\j satisfies 

{I(g)AN)AN{uJi U 0J2) = {I(E)An){An{uji) LJJ An(uj2)) 

= {I(g)AN)AN{uJi) UJ (I(E)AN)ANi0J2) 

= {An<E)I)An{uji) UJ (An(E)I)An{oj2) 
= {An<E)I)An{uji Ut^2)- 

Thus we have estabhshed the structure of a bialgebra. Substituting Defini- 
tion [S] in l|12() yields the recursion for the antipode (|14|l . □ 

Both the definition of the coproduct An and thus the antipode Sn are 
recursive and difficult to use in practice. To develop non-recursive formulae 
for these, it is first necessary to define certain cutting operations. 

Definition 6 For a given forest oj G OF, a parent is any node p with at least 
one branch growing from that node and the children are the nodes branching 
from p. Let Pc denote the number of children of p. Cutting off a child node 
equates to removing the edge connecting the child to its parent. 

— A nodal left cut of degree c is a cut where the c leftmost children of a 
given parent node p are cut off. We can write a nodal left cut as ip{c) 
where < c < Pc- The cut splits a forest uj into two sub-forests, P^i'('^^(w) 
and R^p^'^\uj), where P is the part cut off, with the forest containing the 
c leftmost children of p as root nodes and R is the remaining bottom part 

of LU. 

— A left cut is a collection of < fc nodal left cuts £ — {^pi (ci)}*L]^, where 
{Pi}i=i are distinct nodes of uj. This splits in fc cut-off forests {uji}^^^ 
and a remaining forest R^{lu), where uji is the forest containing the nodes 
connected to the c,; leftmost children of pi and R^{lo) is the forest of the 
nodes connected to the original root nodes. We define P^{lo) ^ N as 

P'^{uj) ^ UJ1UUJ2U ■ ■ -Uujk- (16) 

Note that the definition of a left cut includes the case fc = 0, called the 
empty cut, where R^{uj) — lu and we define P^{lo) — I. 

— An admissible left cut is a left cut, containing the restriction that any 
path from a node in uj to the corresponding root is cut no more than 
once. 

We denote by LC, NLC and ALC the set of all left cuts, nodal left cuts and 
admissible left cuts. To define the coproduct we need to slightly extend the 
definition of an admissible left cut, which we choose to call a full admissible 
left cut. The full admissible left cuts of w e OF are obtained by adding 
an (invisible) root node to form the tree r = B^{uj) G OT, applying an 
admissible left cut on r, and finally removing the invisible root node again. 
We denote by FALC the set of all full admissible left cuts. Thus FALC(w) = 
ALC(t) and for any £ e FALC(w) we have P^(w) = P^(t) and R^{uj) = 
B~{R^{t)). Note that FALC(a;) contains the 'cut everything', where fc = 1, 
P^{lli) = Lo and R^{uj) = I, as well as the empty cut with fc = 0, P^{uj) = I 
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and R^{lo) — uj. It is useful to note that the order in which the cuts are 
performed does not affect P^iuS) or R^{uj). The order of the cuts is taken 
care of by the use of the shuffle product in the definition of P^{uj) given by 
equation Hlt)|) . 

As an example, we list all the cuttings of an example tree in Table |31 
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Table 3 The cuts £i of an example tree r, where £12 is the 'cut everything' full 
cut. Thus NLC(t) = LC(r) = {£0, • • • ,^11}, ALC(r) = {4,..., £5} 

and FALC(t) = {£0, . . . ,^5} U {£12}. 



Proposition 2 The coproduct An o/Hjv is non-recursively defined as 

An{uj)^ J2 P\io)®R^{Lo). (17) 

^eFALC(w) 

Proof To prove that the recursive definition and the non-recursive def- 
inition 117() of the coproduct are identical, an induction argument on the 
number of vertices is used. First recall that YKLG{B~ [t)) — ALC(t), and 
for any i e FALC(S-(t)) we have P^{B-{t)) = P^(r) and R^{B-{t)) = 
B~{R^{t)), this implies that 

^ P\B-{T))®Bt{R^{B-{T)))^ ^ P^{t)®W{t). 

jGFALC(B-(t)) jeALC(r) 
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Using this fact, the coproduct now takes the form 

An{ujt) =UJT(E)1+ i ^ P\u)(g,R^{Lu)\u-i J2 P^{t)®R^{t) 

yeFALC(w) / \jeALC(r) 

= cjT®I+ ^ ^ P\uj)UP'-i{t)®R\u;)W{t) 

£eFALC(w) jGALC(r) 

= ^ P\ujt)®R^{lot). 

^eFALC(wr) 

The last equahty is true because the sum over £ G FALC(ci;) and j e ALC(t) 
is equivalent to the sum over i S FALC(a;T) except for the 'cut everything' 
cut which is equal to the term ujt(E)I. □ 

Corollary 1 The dual of the coproduct Aj^ is o the Grossman-Larson prod- 
uct, that is for any u G N and uji,uj2 G N* we have 

{ijJiOL02,uj) — (a;i®tJ2, ^jv(w)). 

Proof If the sum in H17() had been over ALC instead of FALC, then the dual 
would have been the left grafting. To see this, we use the characterization 
of left grafting in Lemma and observe that the nodal left cut corresponds 
to the dual operation of attaching a number of trees in a given order to a 
common node, while the shuffles in corresponds to the dual operation 

of attaching the forests in all possible ways to different nodes. From Defini- 
tion 13 we see that when the sum is extended from ALC to FALC, then we 
obtain the dual of the Grossman-Larson product. □ 

To present a non-recursive definition of the antipode, we define the rever- 
sal map Sf : N N a.s 

SM^l (18) 
Sf{tiT2 ■ ■ ■ Tj) = (-l)VjT,_i • • • Ti, for all n • • • r,- e OF, 

extended to N by linearity. Thus Sp is the unique anti-automorphism of the 
concatenation algebra which sends r i-^ — r. 

Proposition 3 The antipode Sn ofTi.N is non-recursively defined as 



yeLC(w) / 

Proof In order to prove this result, we need some results about a Hopf alge- 
braic structure of the Free Associative Algebra (FAA) [^1] . Given an alphabet 
A, FAA is the vector space formed by taking all finite linear combinations 
of words over A. In our case, the alphabet is OT, the words are OF and the 
vector space is N . A Hopf algebraic structure Tip \s obtained by taking the 
product ^.p = fi^ as the shuffle product and the coproduct Ap defined as 
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the dual of the concatenation product. The antipode is the map Sp defined 
in 1)18(1 . We need a characterization of Z\i? and Sp in terms of cutting op- 
erations. For a uj G OF let the set of Word Cuts (WC) be a simple cut i 
which splits a word cu into two parts oji = P^{u}) and iu^ = R^{u>) such that 
ix) — ujiuj2- WC contains both the empty cut where P^{uj) = I, R^{uj) = u 
and cut everything where P^{uj) — uj, R^{oj) — I. Note that the difference 
between ALC and FALC is that FALC may contain a nonempty cut from 
WC. A direct definition of zi^^- is 

Apiuj) = P\uj)®R\uj), for all uj G OF. 
^ewc 

From jni we find for w G OF \{I} that 

Q = {SFi^I){oj)= liF{{SF®I)^F{^))= Y Sf{P\uj))UR\uj). 

^ewc(w) 

Thus we find a recursive definition of the antipode Sf 

5f(I) =1, 

Sf{uj) = - Y1 Sf{P\uj))UR\u), (20) 

£eWC(w)\c.o. 

where c. e. denotes cut everything. Now we repeat the same computation for 
S'at, using 1)17(1 . This gives the recursive definition of the antipode Sn 

Sn{uj) = - Y Sn[P\uj))UR\u). (21) 

£eFALC(a))\c.c. 

We prove ((19(1 by induction on the number of nodes. Plugging 1(19(1 into ((21(1 . 
we find for w G OF \{I} that 



Sn{uj) = ^ J2 SfI J2 PHP'i^))uR'iP'{^))]uR'i^) 

<'eFALC(w)\c.o. \jeLC(Pf(w)) / 

= - E E SF{P'{Lo)UP'{R^iu;)))uR'{R^iu;)) 

= Y (^'(^)) ^ f - E Sf{P'{R'{u;))) U R'iR^iu;)) 
jeLC(w) \ iewc{uj}\c.o. 

= E Sf{P'{uj))uSf{R'{oj)) 

= SfI P^H^R^i^) 

\j6LC(a;) 
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We have used the relation Spit^i U 1^2) = »S'_f('^i) LI 6*^(1^2) see Corohary|21 
and the recursion l|20|l , as weU as a careful replacement of the summations 
over £ G FALC(ijj) and j G LC(P^(ijj)) with an equivalent sum over j G LC(w) 
and£GWC(a;). □ 

As an example, we compute An and Sn for the word 1^ = » i ■ The cuts 
LC and FALC are shown in Table0| From the direct formulae we find An{llj) 
and §Ar(ci;) as listed in Table Eland Table El 
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Table 4 Cuts £i of an example word ui. The cuts {£4,£5,ie} are full cuts where 
the leftmost child of the invisible root is cut and £7 is the full cut where both 
the children of the invisible root are cut. Thus NLC(aj) — {^1,^2}, LC(a;) = 
{£o,...,£3}, ALC(uj) = {£o,...,£2}, FALCH = {£0, . . . ,£2} U {£4, . . . ,£7} and 

wcH = {^0,^4,^7}. 



We complete this section by listing some well known but very useful 
relations of Hopf algebras, see Sweedler 37 for further details. 

Corollary 2 '371 Given Hn a Hopf algebra, with product U, coproduct Aj^ 
and antipode Sj\j, then for all Wi, 0^2 G OF 

{Sn<E)Sn)An{uji) = An{Sn{uji)), 

SNitOl) U Sn{uJ2) = SN{i^l LI W2). 

IfTi-N is either commutative or cocommutative, then Sn{Sn{ll})) = uj for all 
uj G OF. 



4 Hopf algebras related to Tijv 

There are two interesting commutative graded Hopf sub-algebras of Hat ob- 
tained by restricting from the set of ordered rooted trees to either the set of 
tall trees (that is trees where each parent has one child) or bushy trees (that 
is trees where there is only one parent). These Hopf sub-algebras are useful 
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respectively for determining the order conditions for the problem iQ), when 
f{y) is constant, or when the numerical scheme has high stage order. 

In this section we will also show that the Hopf algebra He of Butcher, 
based on unordered trees, can be identified as a sub-algebra of Hn- Finally 
we find that the Hopf algebra Hp oi the Free Associative Algebra is related 
to Hat through the operation of freezing vector fields, which can be defined 
as a quotient construction on Hn- 

4.1 Connections to the Butcher theory 

Let T denote all unordered trees and F denote all unordered forests, defined 
as the set of all empty or non-empty unordered words over the alphabet T. 
Recall from JHI the following definition. 

Definition 7 Given the real vector space C — M(F), denote the Hopf algebra 
of unordered forests as He = {C, ^ejUe, ^e,ee, Se)- The product fie ■ 
C®C ^ C is defined as the (commutative) concatenation 

The unit element uc : R ^ C, is given by mc(1) — I- The coproduct Ae ■ 
C C®C is defined by linearity and for any r — B^{uj) E T and u},uj eF 
by the recursion 

Ae{l) =101, 

Aeir) = T®I + {mB+)Ae{B-{T)), (22) 
AeicJT) = Aeicj)AeiT). 

The counit ec : C* — > R is defined by bat (I) = 1 and €^{^1) — ior ui E 
OF\{I}. The antipode Se : C C is, as usual, the two-sided inverse of the 
convolution in He, see TIP for details. 

The main tool used to provide the relationship between the Hopf alge- 
bras of ordered and unordered forests is the symmctrization operator defined 
below. 

Definition 8 The symmctrization operator H : N N is a. map defined by 
linearity and the relations 

= I, 

The shuffle product permutes the trees in a forest in all possible ways, 
and the symmctrization of a tree is a recursive splitting in sums over all 
permutations of the branches. The symmctrization defines an equivalence 
relation on OF, that is, if 
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Thus uji ~ ^2 if and only if ^2 can be obtained from lui by permuting the 
order of the trees in the forest and the order of the branches of the trees. We 
see that an alternative characterization of fl is 

f2{uj) = (t{uj) UJ. 

The integer cr(w) is the classical symmetry coefficient, defined for trees and 
forests as 

a{I) = 1, 

(j{tiT2 • ■ - Tfe) = Cr(ri) • ■ •Cr(Tfe)/il!^2! • • ■ , 

where the integers fii, fi2, ■ ■ ■ count the number of equivalent trees among 
Ti, . . . , Tfe. In other words, if we consider the full group of all possible permu- 
tations of trees and branches acting on a forest uj G OF, then 17(0;) is the size 
of the isotropy subgroup i.e. the number of permutations leaving w invariant. 
The total number of permutations acting on a given forest uj is given by the 
integer n(uj) defined as 

^(I) = 1, 
7r(TiT2 • • • Tfc) = fc!a-(ri) • • • cr(Tfe), 

7r(B+(Tl • • • Tfc)) = 7r(Tl • • • Tfe). 

Note that once a forest lo is symmetrized, then another application of the 
symmetrization yields the scaling 

^tt{u;)Q{lj). (23) 

Let F be the unordered forests. Clearly, there is a 1-1 correspondence 
between unordered forests and equivalence classes of ordered forests, thus 
there is a natural isomorphism F ~ OF / ^. Through this identification, we 
can interpret SI as an injection fl : C ^ N where C = M(F) and N — R(OF). 
From (|23|l we see that the map Q^^ : N ~> C defined as 

^ — ' nlLU I 
wGOF ^ ' 

where forget : ^ C is the natural identification of an ordered forest with 
the corresponding unordered forest, defines a left-sided inverse f2^^{n{b)) = 
b for all beC. 

Theorem 2 The symmetrization operator Q : C N defines an injective 
Hopf algebra homomorphism from the Hopf algebra Tic of unordered forests 
into the Hopf algebra Hn of ordered forests. 
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Proof A Hopf algebra homomorphism is a bialgebra homomorphism, which 
is a hnear map that is both an algebra and a coalgebra homomorphism. f2 
is an algebra homomorphism if 

niucii)) = UNii), 

These conditions are automatically satisfied by Definition|H| Q \s a. coalgebra 
homomorphism if 

eN{n{u)) = ecH, (24) 
AN{n{Lu)) = {Q(g)f2)Ac{uj). (25) 

The first condition 124() follows immediately given that the counits are only 
non-zero when the argument is the empty forest. 

The second relation follows using an induction argument. First we need to 
establish a useful relationship between f] and Ac- Using the sumless Sweedler 
notation 

we find: 

{n(»Q)Ac{uj) LU {n(»Q)Ac{T) = {n®Q) (w(i)®tj(2)) LU (r2®J7) (t(i)®T(2)) 

= l2(tJ(l))®l2(W(2))LiJl2(T(i))®l2(T(2)) 

= ^2(W(1)) U f2(T(i))®l2(tJ(2)) U Q{t(2)) 

= ^?(W(l)T(i))«)/?([J(2)T(2)) 

= (J7(g)/2)w(i)T(i)®u;(2)T(2) 

= {n®Q)Ac{iOT). (26) 

Now we prove H25II by induction. For a forest we find using (|f 5|l and (|26|l 

AN{f2{ujT)) = ziAr(r2(w) U Q{t)) 

= Z\Ar(r2(w)) LU Z\Ar(l2(a;))) 

= {Q®Q)Ac{uj)yiA{n®f2)Ac{T)) 

= (r2®f2)(Z\c(w)Z\c(T)) 

= (r2®l2)Z\c(wT). 

For a tree r = B^{lo) G T and G F wc find using the definition of Q and 
the recursion formulas and H22|) that 

Z\Ar(l2(T)) - Z\Ar(/?(B+(cc>))) = Z\Ar (B+ (r2(^))) 

= B+{Q{uj))®l+{mB+)AN{n{uj)) 
= Q{t)®1 + (I«)B+)(/?cg)l2)zic(w) 
= {n®Q){T®l+ {l®B+)Ac{uj)) 

= {mn)Ac{T). 

The proof extends to a general a G iV by linearity. □ 
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A consequence of the above theorem is a connection between the an- 
tipodes of the corresponding Hopf algebras. 

Corollary 3 Given f2 : C ~* N is a bialgebra homomorphism then for uj G 
F, it follows that 

SNmuj)) = n{Sc{^)). 

Note that symmetrization operator J7 is invertible, so expressions for the 
product, coproduct and antipode of the Hopf algebra of unordered forests 
can be directly expressed in terms of the corresponding functions in the 
Hopf algebra of ordered forests, they are 

Ac{oj) = {n®Q)-^AN{n{uj)), 

Sc{io)^n-\SNmu;))). 

In the final part of this section we will elaborate on the connections be- 
tween the LS-series built on ordered trees, and their commutative counterpart 
the S-series. These series belong respectively to the dual spaces N* and C*, 
are naturally associated through the dual map fi* : N* C* taking the se- 
ries of ordered forests to unordered forests. If a G N* and f3 — f2*{a) G C*, 
we find 

= {n*{a),uj) = {a,n{uj)) = <j{uj) ^ 

On a manifold with a commutative Lie group action the elementary differ- 
ential operators J^{uj) do not depend on the ordering. Thus we find that the 
S-series of /3 as defined in POj equals the LS-series of a as given in 

LS(a) = h^Mu^mc^) - E ^^^^^(^) = m- 

This shows that the normalization Xjaiu}) in the commutative case is com- 
patible with our normalization in the non-commutative case. 

It is interesting to characterize the image of the logarithmic and exponen- 
tial series under fl* . If a is logarithmic (Lemma ISJ then 

m - a(I) = 0, 
/3(wr) = (a, J7(wr)} = (a, n{uj) U I2(t)) = 0, for tJ, r ^ I, 

thus (5 is non-zero only on trees. If a is exponential then 

/3(ti • • ■ Tfc) = (a, I2(ti • ■ • r,-)) = (a, /?(ti) U ■ • • U l2(Tfc)} 
= a(r2(ri)) ■ • • «(l2(Tfc)) = /3(ri) • • • /3(rfc). 



This is a well known condition in the composition of B-series, see [H]. 
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4.2 The Free Associative Algebra and frozen vector fields 

In the proof of Proposition we defined the Hopf algebraic structure Tij? as 
the Free Associative Algebra built of words over an alphabet A, where 
is the shuffle and Ap the dual of the concatenation product. We will briefly 
comment upon the connection between Tijv and Tip in the context of Lie 
group integrators. 

In the theory of Lie group integrators, it is common to call constant 
sections {A4 — > fl) frozen vector fields. If / is frozen, then the Lie derivative 
g[f] — Q for all vector fields g. On the algebraic side, a tree r e OT represents 
a frozen vector field if the left grafting of anything non-constant to the tree 
is zero, 

t[t]=0, for all? G 0T\{I}. 

In this case we see that the Grossman-Larson product becomes just the 
concatenation product tot = tt. The freezing of certain vector fields can 
be understood as the quotient Hp = Hn/G, where G is the linear span of 
any r G OT\{I} grafted to a frozen vector field. As a special example, we 
consider the case where all single-node trees are frozen, so that taller trees 
cannot be produced. Letting A = {B^ {l)}i(zx — {oijigi be the alphabet of 
all single-node trees, we find from H13|l the following well known recursion for 
the coproduct 

Ap{l)= I® I, 
Ap{u)ai) = woi®! + Z\F(w)(8)(I(8)ai). 

5 Concluding remarks 

We have in this paper investigated the algebraic structure of the Hopf al- 
gebra underlying numerical Lie group integrators. We have developed both 
recursive and direct formulae for the coproduct and the antipode, and we 
have in particular emphasized the connection to the Hopf algebra of classi- 
cal Butcher theory and to the free associative Hopf algebra. We believe that 
this work is of particular interest for the construction of symbolic software 
packages dealing with computations involving algebras of non-commutative 
derivations. The algebraic structure of Ti^r is of a universal nature and should 
be of interest also outside the field of numerical integration, for example, in 
the renormalization of quantum field theory and the Chen-Fliess theory for 
optimal control. 
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Table 5 The coproduct in Tijv for all forest up to and including order four. This 
is the dual of the product in the Grossman-Larson Hopf algebra. 
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Table 6 The antipode in Ti.N for all forests up to and including order four. This 
is the dual of the Grossman-Larson antipode. 



